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Abstract 

We consider a directed percolation process at its critical point. The 
probability that the deviation of the global order parameter with re- 
spect to its average has not changed its sign between and t decays 
with t as a power law. In space dimensions d > 4 the global persistence 
exponent 8 p that characterizes this decay is 8 P = 2 while for d < 4 its 
value is increased to first order in e = 4 — d. Combining a method 
developed by Majumdar and Sire JlJ with renormalization group tech- 
niques we compute the correction to 8 P to first order in e. The global 
persistence exponent is found to be a new and independent exponent. 
We finally compare our results with existing simulations. 
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1 Motivations 



1.1 Directed Percolation 

At the initial time A particles are placed randomly with density po on the 
sites of a <i-dimensional hypercubic lattice. They perform independent sim- 
ple random walks with a diffusion contant A. Multiple occupancy is al- 
lowed. The A particles undergo three reaction processes: coagulation upon 
encounter at a rate k, branching at a rate k' , spontaneous death at a rate 7, 

A + A^A 

A^A + A (LI) 



As the branching rate k' is decreased below a threshold value k' c (equal to 7 
in mean-field), the steady state of this system exhibits a continuous transi- 
tion from a state in which a finite positive density of A's survive indefinitely 
to an absorbing A- free state. The order parameter of the transition is pa, 
the average of the local density of A's. 



We have used here the language of the Schlogl reaction-diffusion process 
to describe directed percolation as in Various alternative formulations 
exist (|^, Q). Furthermore the scope of directed percolation reaches far be- 
yond chemical kinetics, as an overwhelmingly large class of nonequilibrium 
systems possessing a phase transition in their steady state fall in the same 
universality class (cellular automata, surface growth, reaction-diffusion pro- 



cesses). This makes of the process Eq. (1.1) a paradigm for non-equilibrium 
systems with a transition in their steady state. Our knowledge of the be- 
havior of the system in the steady state and during the relaxation stages 
rests on numerical simulations (in low space dimensions, d = 1,2) and on 
analytical techniques (short time series expansions in d = 1,2, renormal- 
ization group in d = 4 — e). The critical regime is characterized by a set 
of three independent exponents: the dynamical exponent z, the anomalous 
dimension of the order parameter 77 and the correlation length exponent v. 
Scaling laws for pa can be extracted from special cases of 

PA (t,k' -k' c , Po )=b-^ 1 ^(b 1 ^\k' -k' c \,b-%b^p Q ) (1.2) 

which holds in the limit b — > 00 with the arguments of T fixed. Similar 
scaling relations exist for correlation functions. 



1.2 Global Persistence 

In this article we want to focus on a property that cannot be deduced from 
the knowledge of the scaling properties of correlation functions alone. We 
first define the deviation of the global time-dependent order parameter with 
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respect to its average: 



(t) = lim L~ d / 2 k(x, 

xGL d 



t) 



(riA(x,t)} 



(1.3) 



In Eq. (E3) we denote by n,4(x, t) the number of ^4 particles at site x at time 
t, in a particular realization of the reaction-diffusion process. The brackets 
(...) denote an average with respect to the set of microscopic realizations 
consistent with the initial conditions, the rules Eq. flLlQ and the diffusion. 



We define the global persistence probability as the probability that ^ 
remain of constant sign between and t. Similar quantities have been con- 
sidered (see 0, H) in critical dynamics of magnetic systems, simply being 
the total magnetization. There it was shown that, following a quench from 
a high temperature disordered state to the critical point, the global per- 
sistence probability decays with time as a power law characterized by a 
universal exponent 6 p . In critical dynamics the persistence probability is 
a quantity that appears naturally in the description of the system while it 
relaxes to its equilibrium state. Our motivation for the present work lies in 
the lack of both qualitative and analytical picture of the onset of longe range 
correlations in non-equilibrium systems relaxing to their steady state. We 
believe that the knowledge of the global persistence probability will shape 
our picture of the way the system organizes at criticality. 

The remainder of the article is divided as follows. We first recall in Sec. 2 
the well known correspondence between directed percolation and field the- 
ory. Following Majumdar and Sire Q, it is possible to obtain the global 
persistence probability from a careful analysis of the autocorrelation of the 
global order parameter. This analysis is performed in great detail in Sees. 3 
and 4. In Sec. 5 we turn to the explicit calculation of the persistence expo- 
nent. In our conclusion we compare our results with existing simulations. 



2 Field theoretic formulation 

There are several ways of mapping directed percolation onto a field theory 
(|§,§)- The resulting field theory involves a field ip whose average is the local 
density of A individuals, and a conjugate field if) ; dropping terms irrelevant 
in the vicinity of the upper critical dimension d c = 4 the corresponding 
action reads 



S r. t - / d d x dt [i>(d t + A(<7 - A))V + yV>V# ~i>) ~ P^{t)4> 



(2.1) 



The parameter g can be expressed in terms of the original reaction rates 
fc, 7 and the mass in the propagator is Act = 7 — k' . The action Eq. (|2.1| ) 
is the starting point of the subsequent analysis. Renormalization group tech- 
niques allow us to focus on scaling laws close to or at criticality, during the 
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relaxation process or in the steady state. From here on, as we shall eventu- 
ally focus on phenomena taking place at criticality, we set a = 0. We now 
summarize a few well-known results on the renormalization of the action 
Eq. ( |2.1[) that can be found e.g. in ||. 



One first defines renormalized parameters and fields as follows 

a 2 

i> = VZtl; R , i, = y/Zfa, A = Z~ X Z X \ K , — = Z^ 2 Z- l Z u u^ £ (2.2) 

where p is a momentum scale. From the one loop expression of the two and 
three-point vertex functions one deduces the values of the Z-factors using 
dimensional regularization and the minimal subtraction scheme. They read 

Z = l + -, Z A = l + f, Z u = l + - (2.3) 
The ^-function has the one-loop expression 

0u = ^ = u(-e + 2 7A + 7 - 7u ) = u(-e + 6u) (2.4) 

where we have introduced the Wilson functions 7$ = n d , i = 0, A, u. 
The /3-function has a stable nontrivial fixed point u* = I + 0(e 2 ). Critical 
exponents are then obtained from linear combinations of the ^(u*), e.g. 
z = 2 — 7^ + 7* and rj = 7*. 

We find it convenient to shift ip by its mean-field expression 

M) = — ^— (2-5) 

Therefore the action expressed in terms of the fields <fi = ip — ip m { and cj) = xp 
reads 

s[^] = [\m+K^£—-m Xgp0 2 
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1 + ^p t 2(1 + %p t) 

(2.6) 



We have thus gotten rid of the initial term localized at t = 0. We will use the 
following notation: G^ n,m - ) denotes the (n + m)-point correlation function 
involving n fields 0's and m fields 0's, as defined in Eq. (fO|), and W^ n ' m ^ 
denotes its connected counterpart. The basic ingredients for a perturbative 
expansion are the free propagator G and the free correlator C, defined by 
the zero-loop expression of G^ 1,1 ) and G^ 2 '°\ respectively. We shall need the 
large time behavior of G and C: 

G(k;t',t) = Q(t'-t)^) 2 exp[-X(k 2 + a)(t'-t)] (2.7) 
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C(k;,MO 



t 2 t 2 

1 > Z < 



-Ak 2 (t<+t>). 



(2Ak 2 ) z 



6 (l_e 2Ak2 *<) 

+ 6(2Ak 2 t<)e 2Ak2i < 
- 3(2Ak 2 t<) 2 e 2Ak2t 
+ (2Ak 2 t<) 3 e 2Ak2 '< 



(2. 



in which we have set t < = min{i,t'} and t > = m&x{t,t'}. Note that the 
dependence on the initial density po has disappeared, this is because we 
are focussing on times large with respect to the time scale set by po- We 
have now in hands the building blooks for a perturbation expansion of the 
expectation values of time-dependent observables. 



3 Autocorrelation function 

Our aim is to find the one-loop correction to the function C(t, t') defined by 

C(t,t') = Jd d r W (2 '°\r,t;0,t') = W {2 ' 0) {k = 0;t,t r ) (3.1) 

which merely is the autocorrelation function of the field ip. In order to 
determine C(t,t') we carry out a perturbation expansion in powers of the 
coupling constant g. The first term of this expansion is of course C(k = 
0,t, if). The first non trivial corrections come in six pieces, each depicted 
by a one loop connected Feynman diagram shown in Fig. 1. The explicit 
calculation of these diagrams combined with Eqs. (|2.2|) and (|2.3| ) allows one 
to determine the renormalized autocorrelation function 

Cr(M') = Z- x C(tj) (3.2) 

Using again the shorthand notations t < = min{t,t'}, i> = max{i,t'}, we 
find 



C R (t,t') = \(Xt < ^(^f~ i A 



l--F(t > /t < )+0(e 2 ) 
b 



(3.3) 



which holds for t < ,t > large with t > /t < finite. In Eq. (|3.3| ) the amplitude A 
reads 



.4 



1 + K 



e /1457 vr 2 96, \ nl 2 , 

-( + In 2 + O e 2 

6 V 600 20 25 / y ' 



(3.4) 



and the function F has the expression 

.2 

F(x) 



19 11 tH 96 12 

1 1 ln2H x ■ 

6 10 25 25 



+ ln(l 



x 



101 

~25~ 
23 



5x + x + x 



50 
26 

25 : 



77 3 
— x 6 
50 

1 5 1 

— x° + —— 
25 2hx 



+ ln(l + x 1 
+ 21n 2 x + 31nxln(l - x' 1 ) 



o q 23 a 
- + x + x +X +-X 



25 



1 

25^ 



(3.5) 



+ 4Li 2 (l - x) + Li 2 (l 



^ - l)Li 2 (x- 2 ) - h^x- 1 ) 
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where hi2(x) = — Jq dtln(l — t)/t is the dilogarithm function. The limiting 
behavior of F is found to be 

F(oo ) = 8329 2 2 96 
v ; 1200 5 25 (3.6) 
F(x) X = X -2{x - 1) ln(x - 1) + 0(x - 1) 

We have listed in the appendix the individual contributions to C(t, t') arising 
from the corresponding Feynman diagrams. 



4 Short time expansion 

The result of the previous section Eq. ( |3.3|) for C^(t',i) holds for all times 
t and t' , with t/t' finite, but the limit t <C f is singular. In this section we 
show that for t -C t' the autocorrelation function Cji(t', t) displays power 
law behaviour with respect to both time arguments and determine the cor- 
responding exponents. In this limit the random variable becomes a 
Markovian process for which the persistence exponent may be expressed in 
terms of well-known critical exponents. In the case of the Ising model the 
Markovian approximation for 9 p is already close to the values obtained by 
simulations R, [5| . 

An appropriate method to study the correlation function for t <C t' is the 
short time expansion (STE) of the field ip{r, t) in terms of operators located 
at the 'time surface' t = 0. Since the Gaussian propagator and correlator 
are of the order t 2 for t — > we expect that the leading term in the STE is 
the second time derivative ip of if), i.e. 

#-,i)-Mr,t)> =c(t)^(r,0) + ... (4.1) 

(For t = the second time derivative of the response field ip is equivalent 
to ip.) The function c(t) is a power of t which can be obtained from the 
difference of the scaling dimensions of ip and ip. Naively, c(t) ~ t 2 . 

To compute the scaling dimension of ip in an e-expansion one could deter- 
mine the additional renormalization that is necessary to render correlation 
functions with ip insertions finite. Fortunately, it is possible to express this 
dimension to every order in e in terms of other critical exponents. For the 
initial density pq = oo there is a similarity between directed percolation and 
the semi-infinite Ising model at the normal transition (i.e., for infinite sur- 
face magnetization). In the latter case the short distance expansion of the 
order parameter field near the surface is governed by the stress tensor |8], ^] . 
Due to the translational invariance of the bulk Hamiltonian the stress tensor 
requires no renormalization. Here we look for an initial field which remains 
unrenormalized as a consequence of the translational invariance (with re- 
spect to time) of the stationary state. 

Our argument applies to any dynamic field theory defined by a dynamic 
functional of the form 

/■oo r- 

S[iP, i>] = ^ dtj d d r (MV - Tty, $]) (4.2) 
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We assume that ip satisfies the sharp initial condition ^(r, 0) = pq. For 
directed percolation we have 

r[^,^] = -A^(a-A)V-^#(^-V') (4.3) 
Correlation functions may be written in the form 

/m n 
vy,, $] n j>(fi,u) n ^j'tj) ex p 

i=l 5=1 (4.4) 

where the functional integral runs over all histories {ip, tp} which satisfy the 
initial condition. 

We now introduce a new time variable t — > t' = t+a(t) (with a(t) > — 1 to 
maintain the time order) and the transformed fields ■0' and ip' with ijj'(t) = 
V>(i') and ip'{t) = ip(t'). At lowest order in a(t) the dynamic functional 
becomes 

S[ip, $ = S[ip', ft] - JdtJ d d r a(t)T[ij/, ft] (4.5) 

where a(0) = has been assumed. 

Performing the time shift in the correlation function G^ n ' m ^ and compar- 
ing the terms of first order in a(t) on both sides of Eq. (4.4) one finds 

K i=l 1 5=1 j ) (4.6) 

id. 



= <^(fi,ti)-...^(r n ,t n ) y di J d d ra(t)T[i,i}) 

Here the angular brackets indicate the average with respect to the weight 
exp(—S[iJ),'ij)]). We may choose 

a(t) = a (l-e~ vt ) (4.7) 

to obtain in the limit v — > oo 

E W + E ^ G^ m )({r,t}) = Wr^) • . . . V(r„,i„) / d d rT + ) 
^=1°^ i=i OT i/ J (4.8) 

where T + denotes the operator T[if),t/j] in the limit t —* + . (We have 
assumed that all time arguments of the correlation function are nonzero.) 

This result shows that T + remains unrenormalized to every order of the 
perturbation theory. Therefore its scaling dimension is given by d(T + ) = 
d + z. At the upper critical dimension d c = 4 we find d(T + ) = d(tp) = 
6. In fact, one can show that T + and V>(0) differ for pg 1 = only by a 
constant prefactor. To see this we express T[tjj, ft] in terms of the shifted 
field (j) = if) — ifimf. Since cj)(t),ip(t) ~ t 2 for i — ► while V ; mf(i) ~ £ only 
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the term — (Xg/2)ip^ ai Tp contributes. Thus T+ ~ ?/>(0) ~ ?/>(0), and the STE 
in Eq. ( f4.1| ) becomes 

4>(r,t)-(i;(r,t)) =c(t)T + (r) + ... (4.9) 

with 

c t t \ ^ t -(d(4>)-d(T+))/z _ r ((d+ V )/2-(d+z))/z / 4 

Combining the STE with the general scaling form of the autocorrelation 
function one obtains 

c a (f ■,*) ~ a-*'* (I) (4.ii) 

which holds for t/tf — ► 0. 

5 Global persistence 

5.1 A detour via the Ornstein-Uhlenbeck process 

Let X(t) be a Gaussian process with the following autocorrelation function 

(X{t)X{t')) = e - e v ) \ T '- T \ (5.1) 

for r, r' large (but arbitrary r — r'). The random variable X is thus a Gaus- 
sian stationary Markov process (of unit variance). It satisfies a Langevin 
equation 

^ = -0? ) X(t) + C(t), <)>0 (5.2) 

where £ is Gaussian white noise: 

(C(t)C(t')>=20(°U(t-tO (5.3) 

Therefore X is an Ornstein-Uhlenbeck process. For such a process the prob- 
ability that X be positive between and r decays exponentially as 

Prob{Vr' G [0,r], X(r') > 0} oc e"^ (5.4) 
These are standard results. 

5.2 Expansion around an Ornstein-Uhlenbeck process 

We now consider a Gaussian stationary process X(t) which has the auto- 
correlation function 

(X(t)X(t')} = e-^V-^l + e/(r' - r) (5.5) 
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with /(O) = and eCl. Then X is not a Markovian process. Majumdar 
and Sire [Q] have shown how to evaluate the probability that X be positive 
between and r to first order in e. They find 



Prob{Vr G [0,r], X(r') > 0} oc e"^ 7 



where 



20 



(0) 



7T JO 



dr- 



;i - exp(-24 0) r)) 3 / 2 



+ 0(e 2 ) 



Hakim [10] has extended this result to 0(e). 



(5.6) 



(5.7) 



5.3 Application to the global order parameter 

At any fixed time t there exists a dynamical correlation length £ ~ i 1 / 2 such 
that the system may be considered as a collection of effectively independent 
blocks of linear size £. Hence \I/ is the sum of (L/£) rf independent degrees 
of freedom, which is a Gaussian variable in the limit L — ► oo. We are now 
in a position to apply the result Eq. ( |5.7|) to the random variable 



X(r) = (e r )/yvartf(e T 
which has the autocorrelation function 



(X(r)X(r')) = 
with, after Eqs. ( p~TlD and flO] ) 



e -<V'-r| F(e |r'-r| 
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Substitution into Eq. (5.7) yields 



0? 



where the integral 



i + +*(«>) 



00 x 3 F(a;) 
(x 4 - l) 3 / 2 



(5. 



(5.9) 



(5.10) 



(5.11) 



(5.12) 
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has the analytic expression 
_ 13 9C 91vr 3vr 2 



200 20 

, r(|) 2 



200 



16 



vr 1 41 
In 2 + — 

4 8 80 



1 1 5 3 

+ 50 3F2( 4' 1 ' 1; 4'2 ;1) 




77vr 23 3, n 

1 In 2 

200 12 4 



13 5 5 



5 7 7 



i m) 2 3 . 

30 l 4' 1 '4'4'4' j T(i) 2 



2tt 3 1 3 5 3 

3^2' 4' '4'2' ^ (5.13) 



i r(i) 2 i 7 

240W #,2( 4 ,1,1; 4' 2;1) 



1 \/2^3 3 3 g 5 
150rYl)2 ^ 2( 4' '2 ! 4'2 ; j 



1 3 3 7 

+ g 3^2(4' !> ^ 2' 4' ^ 

=0.630237... 



113 3 



3T(1) 2 . ■ • . 
8 v^F l 4'2' '4'2' ; 



where C denotes Catalan's constant and 3F2 the hypergeometric function 
of order (3, 2). The final result reads 

p = 0(°)(1 + 0.134e + C(e 2 )) = 2 + 0.059e + 0(e 2 ) (5.14) 



6 Discussion 



6.1 Comparison to existing simulations 



Recently Hinrichsen and Koduvely |11] have performed a numerical study of 
one-dimensional directed percolation in order to determine the asymptotic 
behavior of the global persistence probability. In terms of the variable \& de- 
fined in Eq. (|l,3j), they find the following results. For the probability that ^ 
remain negative between and t they indeed find a power law decay charac- 
terized by a universal exponent 9 P that has the numerical value 9 P = 1.50(2). 
However they find an exponential decay for the probability that remain 
positive between and t. The latter assertion is in contradiction with our 
finding that the global persistence probability decays algebraically irrespec- 
tive of the sign of \£. A plausible interpretation for such an asymmetry 
could be the following. On the one hand the global persistence exponent is 
well-defined in the regime in which the system has lost the memory of the 

d—rj 

initial condition. This regime takes place for times t such that t~5z~ po S> 1. 
On the other hand, ^ is well defined in the limit of infinitely large systems, 
and then it is the sum of a large number of effectively independent contribu- 
tions, which, on a lattice of size L, forces L>(~ t l l z . Hence for numerical 
simulations to yield acceptable results, care must be taken that the double 
limit L z S> t S> Po 2z ^ d ^ is satisfied. Whether the simulations in [|ll]] ful- 
fill these bounds is questionable. Finally, in mapping the random process 



10 



*&(t) to X(t) we have assumed that the time interval under consideration 
contains only times large compared with p^ 1 so that the regime in which 
X is stationary be reached. Strictly speaking, we should have defined the 
persistence probability over a time interval [to,t], with t S> to S> p^ 1 ■ In a 
simulation the choice to = leads to a persistence probability that enters 
the asymptotic regime after times t very large with respect to p^ 1 . This 
may be another problem with [O]. 



6.2 Some speculations 



It is interesting to use the e-expansion to speculate on the numerical value 
of dp in low space dimension. We define the improved value of 9 P , which 
we denote by 9p pec , by the product of the actual value of 9^ deduced from 
numerical simulations, and the 0(e) correction given by Eq. ( |5.7|) . We use 
recent simulations results of one and two-dimensional directed percolation 



carried out by Lauritsen et al. [12]. 



d 


z 




7/ 


4 0) 


/□spec 
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1.581 


1.097 


-0.496 


1.316 


1.8 


1.50(2) 
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1.764 


0.734 


-0.409 


1.567 


2.0 




4-e 


z 12 


2 ^ 16 


£ 

6 


9 5e 
z 24 




2 + 0.059e 


> 4 


2 


2 





2 


2 


2 



(6.1) 



In the first two lines the exponents z, v and r\ are taken from ]12j ; the value 
of 0p°^ was obtained using the hyperscaling relation 9^ = 1 + d/2z, and 
that of 9 P in d = 1 is taken from Of course the column 0p pec gives but 
a qualitative estimate of the true 9 p that is supposedly closer to it than that 
obtained by the crude e-expansion. These predictions certainly have to be 
tested against numerical simulations. 



6.3 Final comments 

We would like to add some comments on table (|6.1|). We find that 9 p > 2 
in d < 4 which says that both the average and the variance of the time 
during which ^ is of constant sign are finite. The average time depends on 
a parameter which has the dimension of time, but there is no time scale left 
in our treatment of the persistence probability. Therefore the average time 
must depend on a microscopic scale or on p^ 1 , which we have treated as a 
microscopic scale. This may also explain why the average time is infinite 
for critical dynamics: For zero initial magnetization a time scale such as 
Pq does not exist. In critical dynamics one could also define a persistence 
exponent for the critical relaxation from an initial state with nonzero lo- 
cal magnetization. In that case the persistence exponent would read in the 
Markovian approximation 9^ = l + d/(2z) > 1 (as for the problem we have 
treated in this article). 
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7 Appendix 



We introduce the notations t < = min{t, t'} and t > = max{t, t'} and the 
ratio r = t > /t < . The one loop contributions to the function C(t,t') are the 
following. 



Fig. 1(a) 
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Fig. 1(b) 
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Fig. l(c)+(d) 
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FIGURES 




r 

e f 

Figure 1: One-loop diagrams involved in the expression of C(t,t'). Note that 
diagrams (a,d,f) are the only one loop contributions to (^>(q = 0, 0)^i>(q = 0, t)). 
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